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We have studied the diffractive dissociation into di-jets of 500 GeV/c pi- 
ons scattering coherently from carbon and platinum targets. Extrapolating 
to asymptotically high energies (where t m i n — ► 0), we hnd that when the 
per-nucleus cross-section for this process is parameterized as a = croA 01 , a 
has values near 1.6, the exact result depending on jet transverse momentum. 
These values are in agreement with those predicted by theoretical calculations 
of color-transparency. 
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Color transparency (CT) is the name given to the prediction that the color fields of QCD 
cancel for physically small singlet systems of quarks and gluons [1]. This color neutrality 
(or color screening) should lead to the suppression of initial and final state interactions of 
small-sized systems formed in hard processes [2], often referred to as point-like configurations 
(PLC's). Observing color transparency requires that a PLC is formed and that the energies 
are high enough so that expansion of the PLC does not occur while traversing the target 
[3-5] (the "frozen" approximation). To demonstrate that the interactions of a PLC with 
the nucleons in a nucleus are suppressed compared to those of ordinary hadrons requires 
identifying observables which depend explicitly on the cross-sections of the PLC's. Mea- 
surements of color transparency are important for clarifying the dynamics of bound states 
in QCD [6,7]. 

Perturbative QCD predicts that a PLC is formed in many two-body hadronic processes 
at very large momentum transfer [2,8]. Experimental studies of such processes have failed to 
produce convincing evidence of color transparency [9-11]. However, the momentum transfer 
may not have been high enough and/or the frozen approximation not valid under these 
experimental conditions. Evidence for color transparency (small hadronic cross-sections) 
has been observed in other types of processes: in the A-dependence of J/i/j photoproduction 
[12], in the Q 2 -dependence of the t-slope of diffractive p° production in muon scattering [13] 
(where Q 2 is the invariant mass of the virtual photon and t denotes the negative square 
of the momentum transfer from the virtual photon to the target proton), in the yield of 
non-diffractive p° production in deep inelastic muon scattering [14], and in the energy and 
flavor dependences of vector meson production in ep scattering at HERA [15]. In this paper, 
we report a direct observation of color transparency in the A-dependence of diffractive di-jet 
production by pions. 

The pion wave function can be expanded in terms of Fock states: 

* = Oi\qq) + P\qqg) + j\qqgg) + ■ ■ ■ ■ (1) 

The first (valence) component is dominant at large Q 2 . The other terms are suppressed by 
powers of 1/Q 2 for each additional parton, according to counting rules [16,17]. When the 
relative velocities of all participating particles are nearly light-like, time dilation lengthens 
the lifetimes of these states and "freezes" the partonic content of the pion "seen" by the 
other particles. Bertsch et al. [18] proposed that when high momentum pions hit a nuclear 
target, the physically small \qq) components will be filtered by the nucleus and materialize 
as diffractive di-jets. In a more recent calculation, based on a generalization of the QCD 
factorization theorem, Frankurt et al. [19] proposed that these small \qq) components can 
scatter coherently from nuclei producing high mass di-jets. When the transverse momentum 
of the individual jets with respect to the beam direction (k t ) is large, the mass of the di-jet 
must also be large. Frankurt et al. show that for k t > 1.5GeV/c the interaction with 
the nucleus is completely coherent and u(\qq)N ^di-jets M) is small. This leads to an 
A 2 dependence of the forward amplitude squared for asymptotically high energies (where 

^min 0). 
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As a good approximation, the integrated per-nucleus cross-section for producing these 
high mass di-jets grows as A 4 / 3 . The forward amplitude squared provides a factor of A 2 [19] 
and the integral of the elastic-scattering form factor, m J exp{—(3B^t) dt, contributes a factor 
^4-2/3 ( on assumption that Ro ~ A 1 / 3 and /3, which depends on nuclear density, is the 
same for all targets). This should be compared with a oc A 2 ! 3 typical of normal pion-nucleus 
interactions (for which the shadowing cross-section is approximately the full 7T-JV inelastic 
cross-section). Using more realistic nuclear form factors and accounting for higher twist 
effects will change the A-dependence modestly, and electromagnetic contributions should 
produce negligible effects [20]. 

In this Letter we report measurements of the ^4-dependence of the diffractive dissociation 
into di-jets of 500 GeV/c pions scattering coherently from carbon and platinum targets using 
data from Fermilab experiment E791. We recorded 2 x 10 10 7r~-nucleus interactions during 
the 1991/92 fixed-target run at Fermilab using an open geometry spectrometer [21] in the 
Tagged Photon Laboratory. The segmented target consisted of one platinum foil and four 
diamond foils separated by gaps of 1.34 to 1.39 cm. Each foil was approximately 0.4% of 
an interaction length thick (0.5 mm for platinum and 1.6 mm for carbon). Six planes of 
silicon microstrip detectors (SMD) and eight proportional wire chambers (PWC) were used 
to track the beam particles. The downstream detector consisted of 17 planes of SMDs for 
vertex detection, 35 drift chamber planes, two PWCs, two magnets for momentum analysis, 
two multi-cell threshold Cerenkov counters (not used in this analysis), electromagnetic and 
hadronic calorimeters for electron identification and for online triggering, and two planes 
of muon scintillators. An interaction trigger required a beam particle and an interaction 
in the target. A very loose transverse energy trigger, based on the energy deposited in the 
calorimeters, and a fast data acquisition system allowed us to collect data at a rate of 30 
Mbytes/s with 50/is/event dead time and to write data to tape at a rate of 10 Mbytes/s. 

Data reconstruction and additional event selection were done using offline parallel pro- 
cessing systems. The data for this analysis are selected with the primary requirement that 
at least 90% of the beam momentum is carried by charged particles. This reduces the effects 
of unobserved neutral particles. In addition, this analysis uses events produced relatively 
early in the experiment, before the performance of the drift chambers degraded in the region 
through which the pion beam passed. The offline selection for this analysis was implemented 
after most of the data had already been filtered for other analyses, and we use data taken 
only when all five targets were in place. In all, about 10% of the experiment's integrated 
data set is used in this analysis. 

The JADE jet-finding algorithm [22] is used to identify two-jet events. The algorithm's 
cut-off parameter (m cut ) was optimized for this analysis using Monte Carlo simulations which 
are described below. For each two-jet event we calculate the transverse momentum of each 
jet with respect to the beam axis (A;*), the transverse momentum of the di-jet system with 
respect to the beam axis (q t ), and the di-jet invariant mass, Mj. The di-jet invariant mass 
is related to the quarks' longitudinal momentum fractions in the pion infinite momentum 
frame (x) by simple kinematics: Mj = k 2 /[x(l — x)]. To assure clean selection of high- 
mass di-jet events, a minimum k t of 1.2 GeV/c is required. The selection of clean di-jet 
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events was verified by testing their relative azimuthal angle which for pure di-jets should 
be near 180°. This angle is required to lie within 20° of this value. The size of a \qq) 
system which produces di-jets with k t > 1.5 GeV/c can be estimated as 1/Q < 0.1 fm where 
Q 2 ~ Mj > Ak 2 ~ 10 GeV 2 /c 2 . The distance that the \qq) system travels before it expands 
appreciably, the coherence length, is given by £ c ~ (2p lab )/(Mj — m 2 ) [3] which is ~ 10 
fm for Mj ~ 5 GeV/c 2 . Therefore, we expect that the di-jet signal events selected in this 
analysis evolve from point-like configurations which will exhibit color transparency. 
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FIG. 1. q 2 distributions of di-jets with k t > 1.25 GeV/c from interactions of 500 GeV/c tt 
with platinum and carbon targets. 

The q 2 distributions of the selected di-jet events are shown in Fig. 1. The peaks at small 
q 2 arising from coherent scattering from nuclei are smeared due to missing neutrals and 
detector resolution, but the integrated coherent signals can be extracted with reasonable 
accuracy. Before detector acceptance and smearing, coherent peaks should be produced 
with dN/dq 2 oc exp(— bq 2 ) with b inversely proportional to the nucleus's radius (2.44 fm for 
carbon and 5.27 fm for platinum [23]). Because theory predicts that the A-dependence varies 
with k t [19], the analysis is carried out in three kt regions: 1.25 GeV/c < k t < 1.5 GeV/c, 
1.5 GeV/c < k t < 2.0 GeV/c, and 2.0 GeV/c < k t < 2.5 GeV/c. Altogether, we find 
about 5000 coherent di-jet events in the carbon data set and about 2800 in the platinum 
data set. 



To determine the relative number of coherent events produced in each target, we fit the 
data of Fig. 1 as sums of q 2 distributions of di-jet events produced coherently from nuclear 
targets, of di-jet events produced coherently from individual nucleons but incoherently with 
respect to the nuclear targets, and of background. The shapes of the di-jet distributions are 
calculated using Monte Carlo simulations. We use the LUND PYTHIA-JETSET package 
[24] to generate di-jet events with masses of 4, 5, and 6 GeV/c 2 . This covers the range of 



k t observed in the data. The quark momentum distribution inside the pions is generated 
using an asymptotic wave function [25,26] which is consistent with the data presented in our 
companion paper [27]. Coherent nuclear events are generated with slopes appropriate to 
carbon and platinum. Coherent nucleon events are generated with q\ slopes appropriate to 
the nucleon radius (0.8 fm), truncated at qf = 0.015 (GeV/c) 2 to account for the nucleon 
binding energy. The generated events are passed through a detector simulation and digitized 
to mimic real events. They are reconstructed and analyzed using the same programs used 
to reconstruct and analyze the real data. 
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FIG. 2. kt distributions of simulated di-jets with 4 (slanted lines), 5 (vertical lines), and 6 GeV/c 2 
(horizontal lines) masses normalized and superimposed on the distribution from data taken with 
a platinum target. 

To determine the relative efficiencies for observing di-jet events produced through co- 
herent diffractive scattering from carbon and platinum nuclei, we use Monte Carlo samples 
generated with di-jet masses of 4, 5, and 6 GeV/c 2 . The proportions are adjusted to repro- 
duce the kt spectrum of the data, as observed in Fig. 2. For each kt range, and for each 
target, we fit the g 2 distribution using the signal shapes from the Monte Carlo simulations 
and assuming the background contribution is linear in q\. Fig. 3 shows the results of the fit 
forl.5GeV/c < k t < 2.0GeV/c. The dashed line shows the coherent nuclear distribution, 
the dotted line the coherent nucleon/incoherent nuclear distribution, and the dash-dotted 
line the residual background. This background represents the components of the data which 
are not simulated well, such as badly identified jets. The background's contribution in the 



region of the coherent distribution is small. These fits provide normalization factors between 
the number of simulated events of each kind and the data. 

We derive the relative numbers of produced di-jet events for each target in each kt bin 
by integrating over the diffractive terms in the fits and accounting for the relative efficien- 
cies as described above. The signals from the carbon and platinum targets in any one k t 
range have slightly different mass distributions, and we correct the relative yields to account 
for this. Using the measured target thicknesses, we determine the ratio of cross-sections 
for coherent production of diffractive di-jets from the carbon and platinum targets (which 
received essentially the same beam flux). Theoretical calculations of color transparency at 
asymptotically high energies predict a = <JoA a with a = 4/3 assuming very simple nuclear 
wave functions. At the energy of this experiment, \t m i n \ > reduces the cross-sections for 
coherent scattering on carbon (platinum) targets to 0.98 (0.93), 0.97 (0.87), and 0.94 (0.76) 
times their asymptotic high-energy values for Mj = 4.2, 5.0, and 6.0 GeV/c 2 (masses rele- 
vant for our k\ bins). We extrapolate our calculations of a to asymptotically high energies 
dividing the yields by these factors. The results for each k t bin are listed in Table I. Using 
more realistic wave functions, the predicted value of the asymptotic value of a is 1.45 for 
carbon and platinum targets. Frankfurt et al. [19] predict some dependence of a on kt as 
well. These values, labelled a(CT), are also listed in Table I. 

We have considered the sources of systematic uncertainty which are listed in Table II. 
The degree to which the simulations represented correctly the effect of not including the 
neutral component of the jets is checked by raising the minimum total momentum of charged 
particles from 450 GeV/c to 470 GeV/c. The difference in the final results of a with and 
without this requirement is taken to be the corresponding systematic uncertainty ("effect 
of neutrals"). The uncertainty due to using discrete masses in the Monte Carlo simulation 
is estimated using the difference between results assuming that all the events in a given 
kt range have one mass or another ("discrete masses"). A third uncertainty is assigned to 
the change in yields due to mass-distribution differences in carbon and platinum. We also 
observe some sensitivity to the fitting range used; the associated differences are taken as a 
fourth systematic uncertainty. The total systematic uncertainty is taken by adding these 
contributions in quadrature, retaining the signs when not symmetric. 

In summary, we have measured the relative cross-sections for diffractive dissociation into 
di-jets of 500 GeV/c pions scattering from carbon and platinum targets. Extrapolating to 
asymptotically high energies (where t min — > 0), we find that when the cross-section is param- 
eterized as a = cr A a , a ~ 1.6. This is consistent with expectations based on calculations 
of color-transparency models and is clearly inconsistent with the a oc A 2 / 3 dependence ob- 
served for inclusive 7r-nucleus scattering. We have measured a in three ranges of kt] because 
the uncertainties are large, the results are consistent with no variation, but also with the 
predicted variation. The clear diffractive structure of the signals and variation of the co- 
herent cross-section with A indicate that we have observed the coherent scattering of \qq) 
point-like configurations predicted by color-transparency. 
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FIG. 3. qt distributions of di-jets with 1.5 < h < 2.0 GeV/c for the platinum and carbon 
targets. The lines are fits of the MC simulations to the data: coherent nuclear dissociation (dotted 
line), coherent nucleon/incoherent nuclear dissociation (dashed line), background (dashed-dotted 
line) and total fit (solid line). 



TABLE I. The exponent in a oc A a , experimental results for coherent dissociation and the 
Color- Transparency (CT) predictions. 



kt bin 
GeV/c 


a 






Aa 


a (CT) 


1.25 - 1.5 


1.64 


±0.05 


±0.04 -0.11 


±0.06 -0.12 


1.25 


1.5 - 2.0 


1.52 


±0.09 


±0.08 


±0.12 


1.45 


2.0 - 2.5 


1.55 


±0.11 


±0.12 


±0.16 


1.60 
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TABLE II. The systematic errors 



kt bin 
GeV/c 


effect of 
neutrals 


discreet 
masses 


different efficiency 
for C and Pt 


fit range 
sensitivity 


total 


1.25 - 1.5 


-0.09 


±0.03 -0.06 


±0.02 


±0.02 


±0.04 -0.11 


1.5 - 2.0 


±0.03 


±0.03 


±0.06 


±0.04 


±0.08 


2.0 - 2.5 


±0.06 


±0.05 


±0.06 


±0.07 


±0.12 
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